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Abstract 

Two little explored aspects of Compton scattering of the CMB in clusters are 
discussed: The statistical properties of the Sunyaev-Zeldovich (S-Z) effect in the 
context of a non-Gaussian density fluctuation field, and the polarization patterns in 
a hydrodynamcially-simulated cluster. We have calculated and compared the power 
spectrum and cluster number counts predicted within the framework of two density 
fields that yield different cluster mass functions at high redshifts. This is done for 
the usual Press & Schechter mass function, which is based on a Gaussian density 
fluctuation field, and for a mass function based on a ^'^-distributed density field. 
We quantify the significant differences in the respective integrated S-Z observables 
in these two models. 

S-Z polarization levels and patterns strongly depend on the non-uniform distribu- 
tions of intracluster gas and on peculiar and internal velocities. We have therefore 
calculated the patterns of two polarization components that are produced when the 
CMB is doubly scattered in a simulated cluster. These are found to be very different 
than the patterns calculated based on spherical clusters with uniform structure and 
simplified gas distribution. 

Key words: clusters of galaxies, cosmic microwave background, large-scale 
structure 



1 Introduction 



The pioneering experimental and observational CMB work of our esteemed 
colleague Francesco Melchiorri has paved the way for generations of students 
and colleagues in Italy and elsewhere, who are continuing the mission he helped 
define and lead for more than 30 years. Francesco worked hard to advance 
experimental and observational capabilities to measure the Sunyaev-Zeldovich 
(S-Z) effect in nearby clusters. His expertise and leadership resulted in the 
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successful development of MITO, the measurement of the effect in the Coma 
cluster (De Petris et al. 2002), and the upgrading of MITO to a bolometer array 
(Lamagna et al. 2005). As noted by Rephaeli (2005) in Francesco's obituary, 
his passing marks the end of the founder's era in the growth of the thriving 
Italian CMB community. 

In light of current and future capabilities of the many upcoming S-Z experi- 
ments, we have calculated the power spectrum and cluster number counts in 
the context of a specific non-Gaussian density fiuctuation field. As we discuss 
in the next section, certain observational results point to the possible need for 

excess power on cluster scales. The enhanced interest in non-Gaussian mod- 
els motivates contrasting the integrated S-Z observables of one such viable 
model with those of the standard Gaussian model whose parameters are well 
determined. 

Future sensitive mapping of the S-Z effect in individual clusters will allow mea- 
suring its polarization level. A basic description of some of the S-Z polarization 
patterns in clusters has already been given by Sazonov & Sunyaev (1999). In 
Section 3 we briefiy summarize some of the features of two of the polariza- 
tion components that are induced when the CMB is doubly scattered in a 
(non-idealized) cluster whose dynamical state and gas properties are deduced 
directly from hydrodynamical simulations. 



2 Integrated S-Z Observables in a Non-Gaussian Model 

It has been suggested (Mathis, Diego, & Silk 2004) that some observational 
results possibly point to an appreciable deviation from a Gaussian probability 
distribution function (PDF) of the primordial density fluctuation field. These 
are the relatively high CMB power at high multipoles measured by the CBI 
(Mason et al. 2003) and ACBAR (Kuo et al. 2004) experiments, the detection 
of structures with high velocity dispersions at redshifts z = 4.1 (Miley et al. 
2004) and z = 2.1 (Kurk et al. 2004), and the apparent slow evolution of the 
X-ray luminosity function (as deduced from catalogs compiled by Vikhlinin et 
al. 1998 and MuUis et al. 2003). 

Analyses of WMAP all-sky maps (Komatsu et al. 2003, McEwen et al. 2006) 
do not rule out an appreciable non-Gaussian PDF tail on clusters scales. En- 
hanced power on the scales of clusters may lead to detectable differences in S-Z 
power levels and cluster number counts as compared with those predicted from 
the standard Gaussian PDF. It is of interest, therefore, to explore the conse- 
quences of a non-Gaussian PDF, especially in light of the large S-Z surveys 
that will be conducted by ground-based telescopes and the Planck satellite. 
These will measure the effect in thousands of clusters and map the anisotropy 
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it induces in the CMB spatial structure. 



The non-Gaussian model we consider here has its origin in isocurvature fluctu- 
ations produced by massive scalar fields that are thought to have been present 
during inflation. The model is characterized by a PDF that has a xL distri- 
bution of m CDM fields which are added in quadratures (Peebles 1997, 1999, 
Koyama, Soda & Taruya 1999). 

We have carried out detailed calculations of the S-Z power spectrum and 
cluster number counts predicted in the standard Gaussian and a non-Gaussian 
model with m — 1,2. Below we briefiy describe the principal aspects of our 
work and its main results; a more comprehensive description can be found in 
Sadeh, Rephaeli & Silk (2006; hereafter, SRS). 



2. 1 Calculations 



The integrated statistical properties of the population of clusters involve the 
basic cosmological and large-scale quantities, and essential properties of in- 
tracluster (IC) gas. For Gaussian models these properties have been quan- 
titatively determined by semi-analytic calculations (e.g.. Kaiser 1981, Co- 
lafrancesco et al. 1994, 1997, Molnar & Birkinshaw 2000, Sadeh & Rephaeli 
2004) and from hydrodynamical simulations (e.g., Springel et al. 2001). 



The Press & Schechter (1974) mass function is 

n(M,.) = -F(„-^^«, (1) 



where // = 5c/ o and is the background density. The critical overdensity for 

spherical collapse (which is only weakly dependent on redshift) is assumed to 
be constant, bc= 1.69, and the density field variance, smoothed over a top-hat 
window function of radius R, is 

CXD 

a\R) = J P{k)W\kR)k^dk, (2) 





where P{k) = Ak^T'^{k). Evolution of the mass variance is given in 

'^^^'^ = 9[nuo)] i + . ' 
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where 



and (Carroll, Press & Turner 1992) 



The function in equation (1) assumes the form 

F(;.) = J|e-4f , (6) 



for a Gaussian density field, and the form 

for a density fluctuation field which is distributed according to the Xm model 
mentioned above. 

The highest contrast with the standard model is obtained for m = 1, the case 
we have mostly explored and discuss here, but in order to demonstrate the 
impact of this parameter on PDF tail, we have also considered the case m = 2 
(for details, see SRS). Note that the higher the value of m, the closer is the 
PDF to a Gaussian, as would be expected from the central hmit theorem. 

The CDM transfer function for the Gaussian model is 

T{k) = ^" [1 + 3-89g + (16.1g)2 + {5A6qf + (6.71g)^]"i/^ (8) 

For the the xL model, we adopt the isocurvature CDM transfer function 



l + 215g+(16g)V(l + 0.5g) ^ ^ 



, (9) 



with q = k/{VLmh^) Mpc^^ (Bardeen et al. 1986). We adopt the (now) standard 
ACDM fiat cosmological model, with fl/^ = 0.7, = 0.3, h = 0.7. In the 
Gaussian case the spectral index was taken to be n = 1 with mass variance 
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normalization dg = 0.9. In the non-Gaussian model n = —1.8, and by requiring 
that the present cluster abundance is the same as calculated in the Gaussian 
model, the value erg = 0.73 was deduced the Xi niodel. 

The basic expression in the calculation of the S-Z angular power spectrum is 
Ce^ j j n{M,z)Ge{M,z)dMdz, (10) 

2 M 



where r is the co-moving radial distance, and Ge is obtained from the angular 
Fourier transform of the profile of the thermal S-Z temperature change at an 
angular distance 9 from the center of a cluster, , AT{9). The function Ge, 
which is proportional to (AT)^, is fully specified in, e.g., Molnar & Birkin- 
shaw (2000). In the limit of non-relativistic electron velocities the relative 
temperature change assumes the simple form 



AT 



(0) 



xcoth 



X 



y{0). 



fill 



For an isothermal cluster with a (3 density profile, the Comptonization param- 
eter is 



y{0) = 5 , tan 



P 



1 - {0/pOc)' 

^J i + (W 



(12) 



where Uq, Tq, are the central electron density, temperature, and Tc, and 6c are 
the core radius and the angle subtended by the core, respectively. The gas 
profile is truncated at an outer radius R = lOvc- 

The mass-temperature scaling relation is of basic importance for relating the 
magnitude of the S-Z effect to the cluster mass, and thereby to the mass 
function. We have used the following form for this relation 



T{M,z) = Ti, 



M 



1O15/1-1M0 



(13) 



normalized such that the gas temperature of a IO^^H'^Mq cluster is 8.5 keV 
(at z — 0). The parameter a is usually taken to be 2/3, in accordance with the- 
oretical predictions based on hydrostatic equilibrium; we use a range of values 
that ±10% around 2/3. The uncertainty in the evolution of the gas temper- 
ature is accounted for by taking the two extreme values of and 1 for the 
parameter ip, signifying either no evolution, or strong evolution, respectively. 
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The electron density is 

ne{M,z)c^no-^{l + z)^ (14) 



with the scahng of the gas mass fraction to ~ 10% {e.g. , Carlstrom, Holder 
& Reese 2002). Since we do not yet know the redshift dependence of /, we 
assume / = 0.1 to be roughly vahd throughout the redshift interval (0 — 6) 
considered here. Finally, the core radius is cahbrated according to the simple 
relation 

/ M V^^ 1 
Tc = 0.15/i-^ Mpc. (15) 



Observations indicate a variance of ~ 20% in the value of the core radius, 
which is assumed here, with a corresponding range in the value of the central 
density. 

The number of clusters with S-Z flux (change) above AF,^ is (e.g. Colafrancesco 
et al. 1997) 

N{> AF^) = j r^^dz j n{M,z)dM. (16) 



For a cluster with mass M at redshift z 

AF. = ^^^^(x)yo / Rs{\7-ie\,'yB)-y{m,M,z)dn, (17) 



where 7^ and 7 denote line of sight (los) directions through the cluster center, 
and relative to this central los, respectively, 

yo = 2^no(M, z)To{M, z)rc{M, z), (18) 



and the spectral dependence (of the thermal component) of the effect is given 
(in the nonrelativistic limit; more on this later in this section) by 
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where x — hu/kT. The profile of the effect is given in 



y{\ie\,M,z) = 



tan 



-1 



P 



1 - jO/pdc) 
\ 1 + {0/9 



(20) 



which is the los integral along a direction that forms an angle 9 with the cluster 
center. Rs{\'y — 'y^l, ctb) denotes the angular response of a detector whose beam 
size is given in terms of cxs. Finally, 

AF„ ^ i^lfifl. (21) 



is the fiux weighted over the detector spectral response E{i/). 

For determining S-Z cluster number counts we adopt the PLANCK/HFI de- 
tection limit of 30 mJy. This fiux limit translates to a lower limit of the integral 
over the mass function, equation (16). The results for the S-Z power spectrum 
and number counts described below refer to these two models: (a) IC gas 
temperature evolves with time according to = 1, and (b) no temperature 
evolution, = 0. Power spectra were calculated &i v — 353 GHz, with a 7.1' 
beam size. Number counts were calculated also v — 143 GHz and v — 545 
GHz. 



2.2 Results 



The augmented tail of the non-Gaussian distribution function (compared with 
that of a Gaussian) leads to earlier cluster formation (by virtue of the higher 
probability that an overdense region attains the critical density for collapse) 
and enhanced high-mass cluster abundance. Thus, in the non-Gaussian case 
the S-Z power spectrum reaches higher levels, and its peak shifts to higher mul- 
tipoles, reflecting the higher density of distant clusters. This effect is greatly 
reduced in a non-evolving temperature model (in which the 1 -\- z scaling is 
factored out). 

Power spectra calculated over the rcdshift interval [0, 6] for cases (a) and (b) 
are plotted in Fig. 1 and Fig. 2. In these figures the upper and lower panels 
correspond to cases (a) and (b), and the left and right panels are for the 
Gaussian and x{ models, respectively. The figures show the distribution of the 
power spectrum with redshift (Fig. 1) in A2; = 1 and in (logarithmic) mass 
intervals of log AM = 1 (Fig. 2), respectively. The maximum power levels 
attained in case (a) are ~ 3 • 10~^^ and ~ 6 • 10~^^ for the Gaussian and 
non-Gaussian models, respectively. Both spectra peak at multipoles higher 
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Fig. 1. Power spectrum of the S-Z effect. Shown are levels of power contributed by 
clusters in the redshift range of < z < 6, and at 6 redshift intervals. Left and 
right-hand panels correspond to the Gaussian and Xi models, respectively. Upper 
and lower panels relate to cases (a) and (b), respectively 



than i > 10, 000; in the Gaussian model this is due to the high gas densities of 
distant clusters, whereas in the non-Gaussian model this behavior is attributed 
to the combination of high gas densities and the long high-mass tail. 

Distant clusters are cooler in case (b) than in case (a), due to the redshift 
independence of their temperatures, so that power levels are lower (~ 8 ■ 10"^'^ 
and 5 ■ 10~^^) in the Gaussian and models, respectively. In the Gaussian 
model the power peaks at ^ ~ 6000, reflecting the lower contribution from 
cooler, more distant clusters. In the non-Gaussian model, a relatively high 
abundance of massive clusters at high redshifts leads to sustained high power 
levels on these scales. Note that irrespective of whether or not the temperature 
changes with redshift, it still increases with increasing mass, so even when 
= a more abundant population of hot clusters is present in the model. 
In the model the contribution of relatively high z clusters to the power 
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Fig. 2. As in the previous figure, but here levels of the power spectra are calculated 
for four mass intervals spanning the range W^^h~^MQ < M < 10^^h~^ Mq. 

is quite evident; in case (a) the contribution of clusters lying in the redshift 
range < z < 1 peaks at £ ~ 3000, but the total power continues to rise due 
mostly to clusters in the range 1 < z < 2. The same behavior, though less 
pronounced, is evident in case (b). 

The distribution of S-Z power in four mass intervals is plotted in Fig. 2. Of 
course, the highest mass interval (lO^^/i-^M© < M < lO^^h'^MQ is most 
significant in distinguishing between the signatures of Gaussian and non- 
Gaussian models. Clusters with masses in this range contribute negligibly at 
all multipoles in the Gaussian model (for which most of the power is produced 
by clusters in the mass range 10^^h~^MQ < M < IO^^H'^Mq), whereas they 
dominate the total power up to £ ~ 1000 and i ~ 2000 in cases (a) and (b), 
respectively. At higher multipoles the 10^^h~^MQ < M < IO^^H'^Mq range 
dominates. In the Gaussian model the contribution of low-mass clusters in the 
range IO^^H'^Mq < M < IO^'^H'^Mq is dominant only at the highest multi- 
poles. To contrast the predictions of the two models, the ratio of total level 
of the S-Z power in the Xi model to that in the Gaussian model for cases (a) 
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Fig. 3. A comparison of S-Z power spectrum levels in the Xi ^^id Gaussian models 
for cases (a) and (b). 




ro° lo' 10* 



Fig. 4. Predicted S-Z power spectra at 31 GHz together with CBI and ACBAR 
datapoints. 

and (b) is shown in Fig. 3. For i = 10'^ — 10"^, this ratio increases from 7 to 20 
in case (a), and from ~ 4 to ~ 7 in case (b). 

Perhaps the strongest motivation for a non-Gaussian PDF is the apparent 
difficulty in interpreting the power measured at high i by CBI (Mason et al. 
2003) and ACBAR (Kuo et al. 2004) as (mostly) S-Z induced anisotropy in a 
Gaussian model. The predicted power spectra in the Gaussian and Xi mod- 
els are shown in Fig. (4) together with the CBI and ACBAR measurements. 
The results of the calculations were appropriately adjusted to 31 GHz, the 
observation frequency of CBI. It seems difficult to reconcile the Gaussian re- 
sults with the data, while the predicted S-Z power in case (a) of the Xi model 
matches reasonably well the observed power 'excess' (over the extrapolated 
power measured by WMAP) at i > 2000 (see also Mathis et al. 2004). 

A more direct manifestation of the enhanced high-mass non-Gaussian PDF 
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is apparent in the higher (cumulative) counts in this model, and a redshift 
distribution that is different than in the Gaussian model. In Fig. 5 we show 
the redshift distribution of cumulative S-Z number counts. The plots are ar- 
ranged by the specific model (as in Fig. 1 and Fig. 2). Note that in all models 
considered here only clusters with masses higher than IO^^H'^Mq generate 
fluxes that exceed the adopted detection limit of 30 mJy. 



Total counts in the Gaussian model are ~ 4000 and ~ 3000 in cases (a) and (b), 
respectively, whereas these are ~ 30, 000 and ~ 10, 000 in the corresponding 
cases of the Xi model. In the Gaussian model the largest contribution to the 
number counts (at z = 0) comes from clusters in the mass range IO^'^H'^Mq < 
M < lO^^h'^MfT), whose numbers are at least an order of magnitude higher 
than those in the mass range IO^^H'-^Mq < M < 10^^h~^MQ. In case (a) 
of the Xi model, total counts are dominated by clusters in the mass range 
10^^h~^MQ < M < IO^^H'^Mq, but differences in the relative contribution 
of higher mass clusters are not as pronounced as in the Gaussian case. On 



case a, Gaussian 




10^^<M<10" 

- - 10"<M<10" 
10^^<M<10" 




case b, Gaussian 



raseb, 



-10"<M<10" 
■ 10"<M<10" 




Fig. 5. Cumulative S-Z counts as function of redshift in four mass intervals. (The 
contribution from the lowest mass interval, lO"*^^ < M < 10^^ Mq , vanishes in all 
the cases considered here. 
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modal a, Gaussian model a, 




Fig. 6. Cumulative S-Z counts as a function of redshift. Number counts are shown 
for the three specified frequency channels of the HFI instrument on the Planck 
satellite. 

the other hand, in case (b) in this model the total counts are dominated 
by high-mass clusters at < ^ < 3.5. The cumulative counts at the lowest 
redshifts are slightly higher (by a factor ~ 1.3) in the high-mass range than the 
corresponding contribution from clusters in the mass range IO^'^H'^Mq < M < 
lO^^h-^MQ. Cluster counts are also more broadly distributed in redshift space 
in the non-Gaussian model due to their higher abundance at high redshifts. 

In Fig. 6 we plot the cumulative number counts for two more of Planck fre- 
quency channels, u = 143 GHz and u = 545 GHz (in addition tou = 353 GHz). 
In the non-relativistic electron velocity limit the (thermal component of the) 
S-Z spectral function is g{x). At these frequencies, g{x) = —4, 3.2, 6.7, respec- 
tively, so our results for u = 353 GHz can be scaled accordingly. Adopting the 
same value of the limiting flux (30 mJy), we find that values of the ratios of 
number counts at z/ = 353 GHz and u = 545 GHz are larger in the Gaussian 
than in the Xi model. This can clearly be attributed to the significantly higher 
population of clusters with flux exceeding the limit in the Xi model. Note that 
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an accurate calculation of the S-Z intensity change necessitates a relativistic 
calculation (Rephacli 1995) which results in a more cumbersome expression 
for the temperature-dependent spectral function (for which there are a few 
analytic approximations; see, e.g., Shimon & Rephaeli 2004, Itoh & Nozawa 
2004). The deviation from the nonrelativistic calculation can amount to few 
tens of percent for typical temperatures in the range 5-10 keV. Since our main 
focus here is a comparison between predictions of the two density fields for 
quantities that are integrated over the cluster population (rather than an ac- 
curate description of the effect in a given cluster) , it suffices to use the much 
simpler function g{x). 

The impact of varying a and Vc on the predicted power spectra and cluster 
number counts in the Gaussian and Xi models was explored by SRS. Here we 
only note that the observed mean ranges of these parameters do not introduce 
appreciable uncertainty to blur the significant differences between the predic- 
tions of the two models. SRS have also investigated the consequences of the 
Xi model on cluster formation times and the two-point correlation function of 
clusters. 

In summary, the main objective of the work of SRS has been to compare pre- 
dictions of S-Z observables for two mass functions which differ mainly at the 
high mass end. A x^-distributed field is characterized by a longer tail of high 
density fluctuations with respect to the Gaussian model. This leads to the for- 
mation of high-mass clusters at higher redshifts. S-Z power spectra, number 
counts (and two-point angular correlation functions) yield significantly differ- 
ent results in the Xi primordial density field as compared with the respective 
quantities in a Gaussian model. 



3 Polarization Patterns in a Simulated Cluster 

Unpolarized incident radiation is linearly polarized by Compton scattering 
when the radiation either has a quadrupole moment, or acquires it during 
the (first stage of the) scattering process. The degree of polarization due to 
scattering is proportional to the product of the quadrupole moment and the 
Compton optical depth, r. Thus, the largest polarization levels in clusters are 
expected to be at the few fiK level; nonetheless, the detection of polarized 
CMB signals at this level is presumably feasible (Bowden et al. 2004). The 
prospects for measurements of polarized S-Z signals motivate investigation of 
their levels and patterns. 

Several polarization components due to scattering of the CMB in clusters were 
described in detail by Sunyaev & Zeldovich (1980) and Sazonov & Sunyaev 
(1999). In these treatments simple models for cluster morphology and IC gas 
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spatial profile were assumed. Since polarization levels and patterns strongly 
depend on the (generally) non-uniform distributions of IC gas and total mass, 
on peculiar and internal velocities, and also on the gas temperature profile, 
a realistic characterization of the various cluster polarization components can 
only be based on the results of hydrodynamical simulations. This is particu- 
larly the case in clusters with a substantial degree of subclustering and merger 
activity. 

We began a study of S-Z polarization in non-idealized clusters based on cosmo- 
logical hydrodynamical simulations of the formation and evolution of clusters. 
Here we discuss first results for the polarization structure of the properties 
of the S-Z effect. In particular, we consider the two polarization components 
associated with the thermal and kinematic effects that are generated when 
the CMB is doubly scattered, i.e. the polarization is oc r^. In the first, the 
initial anisotropy arises from the thermal S-Z effect, whereas in the second the 
initial anisotropy is produced by scattering in a cluster with a velocity compo- 
nent transverse to the los, Vt = cf3f The spatial patterns of these components 
can be readily determined only when the morphology of the cluster and its 
IC gas are isotropic. Polarization patterns arising from scattering off thermal 
electrons are isotropic (Sazonov & Sunyaev 1999) in a spherical cluster, while 
patterns of all kinematic polarization components are always anisotropic due 
to the asymmetry introduced by the direction of the cluster motion. 



3.1 Kinematic and Thermal Polarization Components 



Linear polarization and its orientation are determined by the two Stokes pa- 
rameters 



Q— — — [ n^dl [ sin"^ a cos 2ip I {a, ip)dD, 

j Uedl J sin^ a sin 2ilj I {a, ilj)dn, (22) 



167r 



where dl is a length element along the photon path, and a and ip define the 
relative directions of the incoming and outgoing photons. The average electric 
field describes the polarization plane whose orientation is determined by the 
angle 



^ = ^tan-i^. (23) 
When the incident radiation is expanded in spherical harmonics 
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l,m 



(24) 



it can be seen (from the orthogonality conditions of the spherical harmon- 
ics) that only the quadrupole moment terms proportional to /2,2 and I2-2 
contribute to Q and U. 

The expressions for the Q and U in the equation (22) are given in terms 
of the relative directions of the ingoing and outgoing photons. Due to the 
obvious dependence of the angles and frequencies on the relative directions of 
the electron and photon motions, these expressions are transformed to a frame 
whose Z axis coincides with the direction of the electron velocity, the axis with 
respect to which the incoming and outgoing photon directions are defined (in 
accord with the choice made by Chandrasekhar 1950). Since the Doppler effect 
depends only on the angle between the photon wave vector and the electron 
velocity, the calculation of the temperature anisotropy and polarization is 
easier if we assume that the incident radiation is isotropic in the CMB frame. 
With the polarized cross section written in terms of angles measured with 
respect to the electron velocity (after averaging over the azimuthal angles). 



da 3aT 



l-^il)PM), (25) 



where P2{f^o) is the second Legendre polynomial, the expression for the Q 
parameter of the scattered radiation is 



Qili) = ^r(l - fil) J P2{fi',)I{fi',)dfi', . (26) 
-1 

Here, /Iq — cos^q, /Iq = cosOq, are the angle cosines between the electron 
velocity and the incoming and outgoing photons, respectively. 

The degree of polarization induced by the kinematic S-Z component was de- 
termined by Sunyaev & Zeldovich (1980) in the simple case of uniform gas 
density. 



where (as defined in the previous section) x — hu/kT. A more complete cal- 
culation of this and the other polarization components was formulated by 

Sazonov & Sunyaev (1999). Viewed along a direction n = {6, 0), the tempera- 
ture anisotropy at a point {X,Y,Z), AT{X,Y, Z,6,(f)) generates polarization 
upon second scatterings. The Stokes parameters are calculated from equation 
(22), 
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Q{X,Y) = ^JdZn,{X,Y,Z) 

X J dil sin^{e) cos(2(/)) A/(X, Y, Z, 6, ^) 
U{X,Y) = ^J dZne{X,Y,Z) 

X Jdn sin'^ie) sin(20) A7(X, Y, Z, 9, 0) (28) 

where now the Z direction is chosen along the los. The intensity change re- 
sulting from first scatterings is 

Ai{x, Y, z, e, <j>)/h = / di{x', y, z', e, 0) 

xn,(X', Y', Z', e, (t>)n ■ (3(X', Y', Z') (29) 
and the optical depth through the point {X, Y, Z) in the direction {9, 4>) is 

r{X, y, Z, 9, (l))^aTl n,{X' , Y\ Z')dl{X', Y\ Z\ 9, 0). (30) 

Q{X, Y) and U{X, Y) fully describe the 2-D polarization field. 

The anisotropy introduced by single scatterings off thermal electrons is the 
usual thermal component of the S-Z effect with intensity change Alt. Using 
the analytic approximation to the exact rclativistic calculation for Alt that 
was obtained by Shimon & Rephaeli (2004), the polarization can be calculated 
as described above 



AIt{X,Y,Z,9,<P) ^'^^Y.Uix) [ ne{r)Qivydl, (31) 

- 1 ^ J 

where = kTe/mc^ and Tg is the electron temperature. In Equation (31) the 
integration is along the photon path prior to the second scattering, fi{x) are 
spectral functions of x that are specified by Shimon & Rephaeli (2004). 



3.2 Clster Simulations 



We investigated the polarization properties induced by scattering in a cluster 
with (a total) mass of 1.4 x 10"^^ Mq simulated with Enzo, an Eulerian adaptive 
mesh cosmology code (O'Shea et al. 2004). The code solves dark matter N- 
body dynamics with the particle mesh technique, and Euler's hydrodynamic 
equations with the piecewise parabolic method algorithm modified for cluster 
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simulations by Bryan et al. (1995). The simulation was initialized with dark 
matter in a volume of space 256h~^ Mpc on a side at 2; = 30 (when density 
fluctuations are still linear on this scale), using a 256^ root grid with 256^ dark 
matter particles (with mass M^m = 1-19 x 10^^ Mq). The standard ACDM 
cosmological model was used with Q,tot = 1, ^^a = 0.7, Q-m — 0.3, Hq — 70 
km/sec/Mpc, = 0.9, and n = 1. Adaptive mesh refinement was turned on 
with an additional four levels of refinement (doubling the spatial resolution at 
each level) and the simulation was run to 2; = 0. At this point the simulation 
was stopped and the Hop halo-finding algorithm (Eisenstein & Hut 1998) was 
used to find the most massive dark matter halo in the simulation. This is our 
simulated 1.4 x 10^^ M© cluster. 

The simulation was then re-initialized at 2; = 30 with both dark matter and 
baryons {Qi, — 0.04, Qcdm = 0.26) with a 128^ top grid and 2 static nested 
grids covering the volume where the cluster forms, giving an effective root 
grid resolution of 512^ cells {0.5h^^ Mpc comoving). With adaptive mesh re- 
finement a maximum spatial resolution of 15.625 h^^ kpc (comoving) was 
attained. This simulation was then evolved to z = 0.06, following the evolu- 
tion of the dark matter and using adiabatic gas dynamics. The cluster was 
then relocated with the Hop algorithm at 2; = 0.06, 0.25, 0.5 and 1.0, and 
data cubes with 256^ cells containing Pdm, Pgas, T^. and the three baryon ve- 
locity components were extracted at these redshifts. The total cube length is 
4/i~^ Mpc, corresponding to a spatial resolution of 15.625/i~^ kpc (comoving). 
The results discussed here are based on the simulated cluster at 2; = 0.5. The 
baryon density and velocity field (with velocity dispersion of ^ 1000 km/sec) 
are shown in Figure 7, and the cluster projected temperature and optical depth 
are mapped in Figure 8. 



3.3 Polarization Maps 



Specific calculations of the above double polarization components necessitate 
ray tracing of doubly scattered photons. In Cartesian coordinates, equations 
(28) & (31) can be written as 



g(X,y) = ^ / I dZdh'n,{v)n,(v') 



^ ^ ^7— E M^M^y 

1^ '-I i=i 



'^Wl') = ^/ /(iZ<<Vn.(r)n,(r') 

^{X^XWzn±f,^,^e(r'Y (32) 



i=l 
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Vetcit)' fEkJ on th? background of baiyon m^s density 




-6. -+ -a -2 -1 Q 1 2 a 4 5 

arcmin 



Fig. 7. The velocity field is shown on the background of the baryon mass density; 
arrow length is linearly proportional to velocity magnitude). Color density scale is 
in units of g/cm? . 

where the first integration (over Z) is along the los. To carry out the compu- 
tationally intensive 4D integrals over the simulation 256'^ data cells, we write 
the last equation in the form 



= 4 j dZn,{v)h+^^{v) (33) 



where P+ and Px are the Stokes Q Sz U parameters, respectively, and 
/+,x g is a convolution of the functions /+,x with g, where 



/+ = ^4 ' fx = 9 = E ^e(r)/.(x)e(r)\ (34) 

i 

The 3D FFT, which can be handled considerably more efficiently, can now be 
calculated using the latter expressions for the core region of the simulation. 

Polarization maps on the sky plane are shown in Figures 9 & 10, with polar- 
ization levels specified in terms of the equivalent temperature. Unlike the case 
of spherical (Sazonov & Sunyaev 1999) or quasi-spherical (Lavaux et al. 2004) 
clusters, we obtain a more intricate polarization patterns and the quadrupole 
structure, which is absent in the spherically symmetric case, is evident. This 
pattern reflects the non-uniform gas distribution and an appreciable degree 
of non-sphericity of the cluster, which substantially affect primarily the dou- 
ble scattering polarization signals, essentially due to the fact that the second 
scattering is more likely to occur far from the cluster center (see also Sazonov 
& Sunyaev 1999, Lavaux et al. 2004). The level of polarization at the 545 GHz 
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Fig. 8. The cluster temperature (color scale in K) and optical depth maps are shown 
in the left and right figures, respectively. 

Planck frequency band, with a 4.5' FWHM beam, can reach ~ 80 nK, but 
when convolved with a narrow beam 1' FWHM profile, its peak value can be as 
high as ~ 0.7 ^K. This level definitely grazes the threshold of next generation 
experiments. It should be noted that we have considered a rich cluster with 
maximal optical depth of tq ^ 0.01, but since this component is proportional 
to ~ r^, typical values could be a factor ~ 4 lower or higher than the specific 
values quoted here. Also, this component was calculated at 545 GHz; lower 
levels of polarization are expected at lower frequencies. 

Equations (28)- (30) for the double scattering kinematic component (oc r^/3) 
can be written in the same form as equation (33), but with 



3 

h+,x = /+,x -*^9 = Y1 /h,x * 9i (35) 



where 



Ji+ = ^ , fix = — ^ — , gi = n{r)(ji{r) (36) 

with Xi [i = 1,2,3) the components of the vector {X,Y,Z). The results for 
the component (for the central region of the simulation) are summarized 
in Figure 10. The level of polarization in absolute temperature units is ~ 1 — 2 
orders of magnitudes smaller than the oc r^G component. Convolving with 
the Planck HFI and a Gaussian beam profile with FWHM of 1', the peak 
values are 3 nK and 25nK, respectively, comparable to the level of the cosmo- 
logical quadrupole polarization (section 3.1). As previously noted, both these 
components are frequency- independent, a fact that has obvious implications 
on the feasibility of separating these components out when the measurements 
will eventually reach the requisite sensitivity. 
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Fig. 9. The r^G polarization component convolved with the Planck HFI 4.5' FWHM 
beam (left) and with a FWHM 1' beam (right). Color scale is in /uK. 



T^^pnlarlialtantorTporarlfPldrtkHFI) ^ ^(j-l rp comporerrt iRVHM 1 iiicmirbaam} 




Fig. 10. The polarization component convolved with the Planck HFI 4.5' FWHM 
beam (left) and with a FWHM 1' beam (right). Color scale is in ^K. 

We end with a summary of the main result of this section: Most polarization 
signals induced by scattering of the CMB in clusters are substantially smaller 
than 1 fiK. An exception is the polarization induced by the double scattering 
thermal effect. Convolving the predictions for this component in the simulated 
cluster with a FWHM 1' beam, we find that the effective polarization signal 
is about ~ 1 fiK. The much improved CMB observational capabilities - with 
projected sensitivities around 1 fiK and arcminute angular resolution - will 
make it feasible to measure the dominant polarization signals in clusters by 
upcoming experiments, thereby supplementing X-ray and (total intensity) S-Z 
measurements to probe IC gas properties and cluster morphology. 
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